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Abstract. Exact solutions of the Ginzburg-Landau (GL) equation for the straight 
film subjected at its edges to the Robin-type boundary conditions characterized by 
the extrapolation length A are analyzed with the primary emphasis on the interaction 
between the coefficient (3 of the cubic GL term and the de Gennes distance A and 
its influence on the temperature T of the strip. Very substantial role is played also 
by the carrier density Ug that naturally emerges as an integration constant of the GL 
equation. Physical interpretation of the obtained results is based on the ris-dependent 
effective potential Vg/ / (r) created by the nonlinear term and its influence on the lowest 
eigenvalue of the corresponding Schrodinger equation. In particular, for the large 
cubicities, the temperature T becomes A independent linearly decreasing function of 
the growing /3 since in this limit the boundary conditions can not alter very strong 
Veff. It is shown that the temperature increase, which is produced in the linear GL 
regime by the negative de Gennes distance, is wiped out by the growing cubicity. In 
this case, the decreasing T passes through its bulk value at the unique density ni°^ 
only, and the corresponding extrapolation length Kt=Tc is an analytical function of 
/3 whose properties are discussed in detail. For the densities smaller than n§'\ the 
temperature stays above Tc saturating for the large cubicities to the value determined 
by ris and negative A while for > rig^ the superconductivity is destroyed by the 
growing GL nonlinearity at some temperature T > Tc-, which depends on A, and 
(3. It is proved that the concentration n^^^ transforms for the large cubicities into the 
density of the bulk sample. 
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1. Introduction 

For more than sixty years, the predictions based on the analysis of the Ginzburg-Landau 
(GL) equations [1] for the order parameter ^(r), supercurrent density j(r) and magnetic 
field B 

— (-iW-gA)^^ + /3|^|^^ + a^ = (la) 
2m 

j = i^(**V* - W**) - — A|*P (16) 

2m ^ ^ m ' ' ^ ^ 

VxB = /xj, (Ic) 

have been in a very good agreement with the experiment [2-4]. In the above equations, 
m is a mass of a superconducting carrier, q is its negative charge, /x is a magnetic 
permeability, A is a vector potential through which the magnetic field is expressed 
according to 

B = V X A, 

and a and ^ > are GL parameters. Developed several years later microscopic Bardeen- 
Cooper-Schrieffer (BCS) theory [5-7] had shown that the superconducting quantum 
(commonly called after this a Cooper pair) is actually two electrons bounded to each 
other by their interaction via phonons with its total mass being equal to the double 
electron mass rUg., m = 2me, and negative charge of the double elementary charge e, 
q = — 2|e|. In fact, the choice of mass in the above equations is arbitrary [6] while the 
double negative charge reflects a fundamental fact of pairing of the two electrons with 
the opposite spins [6]. It is postulated that the density ^^(r) of the superconducting 
electrons is determined by the order parameter ^'(r) according to 

ns(r) = |vl/(r)|V2. (2) 

The expressions for the GL coefficients can be derived from the microscopic theory of 
superconductivity [7,8]: 

^^-^-^ (3) 



2771^ (T) 



iV(0) 



2m^2 (Q) 



2- (4) 



Here, A^(0) is the density of states at the Fermi energy, ks is Planck constant, T^. is 
a bulk critical temperature of superconductor, and ^(T) is a coherence length varying 
with the temperature T as 



e(r) = ^(0)^/^^. (5) 

It means that the cubicity /3 docs not depend on the temperature while the coefficient a 
linearly grows with it, a ~ {T — Tc). For the 'dirty' metals (alloys), the value of ^^(0) in 
(3) and (4) has to be replaced by with I being the mean free path of the material 
[6,7,9]. 
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For the spatially confined samples, equations (1) have to be supplemented by 
the boundary conditions at the superconductor three dimensional (3D) surface, or 2D 
contour, or ID interval S. For the order parameter, de Gennes proposed to use the 
following requirement [6]: 

1 



n(v->iA)* 



A* 



(6) 

s 

with n being an inward unit vector normal to the interface. Analysis of (16) and (6) 
manifests that the real value of the extrapolation length (or the de Gennes distance, as 
it is customarily called on the superconductor slang) A guarantees that no supercurrent 
flows through the confining surface^: 

nJlim(A)=o = 0- (7) 
Numerical magnitude of A characterizes the intensity of the interaction of the 
superconductor with the ambient environment and, as a consequence of this, the 
penetration of the order parameter into it; in particular, for the border with the 
vacuum or the insulator its approaches infinity, takes finite positive values for the 
contact with the normal nonsuperconducting metal and turns to zero while bordering 
the magnetic material. Utilizing BCS equations, it can be shown that, for example, for 
the superconductor half-space terminated at a; = 0, the fading propagation of the order 
parameter into the neighbouring medium at x < is described as [12, 13] 

^{x)r^j e-l^l/^, a; < 0. (8) 

Note that, according to (8), the Cooper pairs can penetrate neither into the vacuum, 
A = oo, (what is natural) nor into the magnetics, A = 0, as expected, since the 
spins of all charged carriers in the magnetic material are forced to align in the 
same direction while in the Cooper pair the spins of each constituent electron are 
antiparallel. Interestingly, for the border with the other superconductor with higher 
Tc the extrapolation length takes negative values what physically means a surface 
enhancement of the superconductivity. This boundary-induced increase of the critical 
temperature was indeed observed in cold worked Ino.993Bio.007 foils [14] and tin samples 
[15,16]. For the former study it was estimated from the increase of the critical 
temperature that the negative extrapolation length is in the order of micron [17], A ^ — 1 
fim, while for the latter case it was argued that the magnitude of the de Gennes distance 
can be controlled by, for example, manipulating the abrasive grain size and the annealing 
time and temperature [15]. Another method of varying the extrapolation length was 
proposed recently from the calculation of the influence of the external electric field £ on 
superconductors; namely, it was argued [18, 19] that its application perpendicular to the 
interface transforms the zero- field de Gennes distance ^£=0 ^^^^ total counterpart 

1 1 



Atot A 



£=0 



I Resonant features of the supercurrent j crossing the interface S for the nonzero imaginary part of A 
have been predicted recently [10, 11] in the framework of hnearized GL equations. 
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with the voltage Ug being expressed through the parameters of the GL theory: 

Us o m ris mc^ 

Here, dimensionless GL parameter k is the ratio of the zero-temperature London 
penetration length A(0) to the coherence length, k = A(0)/.^(0); is superconductor 
ionic background permittivity, and c is the speed of light. While for the classical (low-T^) 
materials with the small k the magnitude of Us reaches high values of ~ 10^ V, for the 
high- Tc superconductors with the large GL parameter, k ~ 100, it drops significantly to 
the value of several tens of volts [18]. 

In the absence of the magnetic fields, equations (la) and (6) simplify to 

A* + /3|^|^* + a* = (11) 

2m 



\s 



T = T, 



E 



(15) 



(12) 

s 

Remarkably, equation (11) is formally identical to the Schrodinger equation describing 
the motion of the quantum particle with the energy 

E=-a= — ^ fl - -1 (13) 
2me (0) \ TJ ^ > 

in the domain confined by the surface S with the effective potential Ve//(r) of the form 

K//(r)=/3|*(r)^ (14) 

Prom (13) one sees that getting the lowest eigenvalue -Emm of the energy is equivalent 
to finding the temperature T of the sample: 

2meM0)" 

Thus, a solution of the problem of the maximal possible minimization of E^^^ shows the 
way of increasing T. The procedure of determining the temperature through the lowest 
eigenenergy Emin of (lo) and/or (11) is widely used in the superconductor research 
[6,20] and will be implemented below too. Note that the effective potential Ve//(r) is 
always nonnegative and, accordingly, taking into account of the cubic term leads to the 
decrease of the temperature T from (15). This has a natural physical explanation: 
the superconducting carriers density n^, which is proportional to the square order 
parameter too, grows with the decrease of the temperature and, accordingly, the role of 
the nonlinear term becomes more prominent as well. 

Armed with the experimental and theoretical knowledge discovered during last 
fifteen-twenty years, we readdress here the problem of the GL description of the 
properties of the straight superconducting film in the absence of the magnetic fields, 
B = 0, partially discussed in the mid 1960s [13] and in 1993 [21,22]. Even though some 
basic equations for this geometry have been derived before, their analysis in the just 
cited references was limited to the linear case and nonnegative A only. In particular. 
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it was shown that the critical film width dc at which the superconductivity emerges, 
depends on A as [13, 21, 23] 

dc = 2^(T)arctan^-^. (16) 

Linearized GL theory at the negative A was also discussed recently for the same 
configuration [17,24-26]. Below, we provide a complete analysis of the temperature 
T dependence on the cubicity f3, superdensity Ug and the extrapolation length A. For 
achieving this, we analytically solve (11) with the GL coefficient a substituted by the 
energy E, equation (13), and, after imposing on its eigenf unction boundary condition 
(12), arrive at the transcendental equation for the calculation of the eigenvalue Emin, 
which, according to (15), defines the film temperature T. Our primary interest is an 
investigation of the combined influence of the extrapolation length A and cubicity f3 on 
the temperature T. It appears that the crucial role in this case is played also by the 
density Us- Special attention is paid to the case of the negative de Gennes distance when, 
as it is proved below, the increasing cubicity promotes the energy from the negative area 
to the positive values for one specific density Ug only. Critical negative extrapolation 
length Ao at which the energy turns to zero, is calculated analytically as a function of 
(3 and its mathematical and physical analysis is performed. 

The results are presented in the following outhne. In section 2, a necessary 
formalism applied to the study of the straight GL strip is derived along with several of 
its asymptotic cases some of which have been known before. Results of the calculations 
for the positive and negative extrapolation lengths are presented separately in section 3. 
The discussion is wrapped up in section 4 by some concluding remarks. 

2. Model and Formulation 

Consider a straight 3D infinite superconducting film of width d. Origin of the Cartesian 
system of coordinates coincides with the middle of the film and the horizontal x axis is 
perpendicular to its edges. Left (right) confining surface is characterized by the uniform 
along the interface de Gennes length A_ (A+). Frequently used method in the study 
of the superconducting systems implements the scaling where the order parameter ^ is 
expressed in units of a//3/|q;| and all distances are measured in units of the coherence 
length ^(T) = ^n?/{2m\a I ) [6, 7, 21] . Advantage of such an approach hes in the fact that 
no any coefficients appear in the resulting differential GL equation and the eigenvalue 
transcendental relationship does contain the film width d [21]. However, shortcomings 
of this method are a continuation of its strong sides; namely, by eliminating the physical 
parameters from the GL equation, such scaling shadows the explicit interrelations 
between them and at the final stage one needs to return again to the units where 
these quantities re-emerge what, due to the quite complicated structure of the derived 
dependencies, is not an easy task to do [21]. For example, in this approach the parameter 
(3 is dropped out from the integration constants of the GL equation and, accordingly, 
it is not clear how to explicitly obtain the influence of its magnitude on the properties 
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of the superconductor and to track the transformation from the hnear to the cubic 
scenario. Since in the present research we are interested in the temperature dependence 
on, among other parameters, the cubicity /3, we employ the scahng where the distances 
are measured in units of the film width d, which is the only natural (temperature 
independent) length for the geometry under consideration. Accordingly, the energy E 
through which the temperature T is expressed according to (15), will be measured in 
units of the ground state energy n'^h'^/ {2md'^) of the ID Dirichlet well of the width d, 
and cubicity j3 will be expressed in units of d^. In this way, the temperature T (via 
the energy E) and cubicity ^ explicitly stay in the differential equation and enter its 
solution and the resulting transcendental equation what allows to directly investigate 
the correlations between these two parameters. 

In the absence of the magnetic field, B = 0, the motion along each of the axes is 
independent from the movement along the other ones. Mathematically, this separation 
of the transverse and two longitudinal motions in the film is expressed via the following 
ansatz: 

^{x,y,z) = e'^^yy+^^'\{x). (17) 

After substitution of (17) into (11) one arrives at the differential equation for the 
transverse function x(a;): 

X" {x) + {^E -kl- kl) X (x) - Px' {x) = 0. (18) 

The designations of the absolute values in the cubic term have been dropped out from 
(18) since it is known that even in the presence of the magnetic fields a solution of the 
GL equations in the simply connected region can be chosen real [7]. As we are interested 
in the lowest possible value of the energy that, according to (15), guarantees the highest 
attainable temperature T, we put ky and kz equal to zero to obtain: 

x!'{x)+Ti^Ex{x)-(5x'{x)=Q. (19) 

Second-order nonlinear differential equation (19) is supplemented by the two boundary 
conditions at the film edges: 



X' {x) ± -^x {x) 



= 0. (20) 

x=±l/2 

Below, for simplicity, we will use the same extrapolation lengths at the opposite surfaces, 
A+ = A_ = A. (21) 



General solution to (19) reads 



x{x)^jlvCsn{Cx + C,rj) (22) 



with 



C = ((£,») = ,/t^. (24) 
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Here, sn{u,k) is a Jacobi elliptic sine [27,28], the integration factor k is determined by 
(or, vice versa, determines) the density of the Cooper pairs, as it directly follows 
from comparison between (22), (23) and (2), and the coefficient C is found from the 
symmetry properties of the film; namely, for the equal de Gennes lengths, equation 
(21), one has either even, Xev {-x) = Xev (x), or odd, Xodd {-x) = -Xodd (x), solutions. 
RecaUing properties of the elliptic functions [27-29] , one immediately gets: 

Cev =K(77) (25a) 
Codd = 0, (256) 

where K(/c) is a complete elliptic integral of the first kind [27-29]. Accordingly, the 
order parameter becomes 

Xev {x) = VC cd{Cx, r]) (26 a) 




'2 

Xodd i^)^ ]Jp sn(Cx, T)) (266) 

with cd(-u, k) being another Jacobi elliptic function (we follow commonly adopted [27- 
30] Glaisher designations [31, 32]). Applying boundary conditions (20) to the transverse 
function x{x), one arrives at 

C (1 - v') sn(C/2, v)-^ cn(C/2, rj) dn(C/2, rj) = (27a) 
for the even states, and 

Ccn(C/2, rj) dn(C/2, rj) + j sn(C/2, rj) = (276) 

for the odd ones with cn(-u, k) and dn(-u, k) being another two major elliptic functions 
[27-29]. These equations determine energy ii^ as a function of the de Gennes distance A 
and, through the factor 77, of the cubicity /3 and parameter k: E = E (A, /3, k). It is the 
aim of the next chapter to analyze these dependencies for the positive, zero and negative 
extrapolation lengths. Since the lowest energy level determining, according to (15), 
the temperature T of the sample, is always even, below mainly this case is discussed. 
It is important to state that the mathematics involved in the derivation of (26) and 
(27) automatically introduces another important physical parameter; namely, it is the 
integration constant k that plays, as we shall see below, a crucial role in determining 
temperature of the film and a square of which, as pointed out above, represents the 
density rig. Note that the coefficients p and k are not completely independent, as 
it directly follows from (23), (24) and (27); for example, if either of them is zero, 
equations (27) do not depend on the second one. Mathematically, this is explained by 
the dependence of the effective potential Veff{x) from (14) on k via the order parameter 
x{x). This has a clear physical meaning too; namely, for the linear case, /3 = 0, the 
density of the superconducting pairs is vanishingly small [6] what corresponds to the 
coefficient k tending to zero, k — )> 0, and, as a result, it is dropped out from (27). Below 
we will interchangeably operate with the factor rj or the parameters f3 and k keeping in 
mind that, if required, the coefficient k, for example, for each given /3 can be determined 
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from (23) after finding a solution E of equations (27). Utilizing properties of the Jacobi 
elliptic functions [27,29,30], it can be shown that the transformation r] ^ l/rj leaves 
{21b) intact while its even counterpart (27a) changes sign of A what means that it is 
sufficient to consider the unit interval 

< ?7 < 1. (28) 

At its right edge, equation (27a) becomes: 
4 1 

we" = . (29) 

Kl-rf 



wit u — 'KyjEjl and, thus, the energy E in this hmit is expressed via the Lambert W 
function [33] 



2„,o^4 1 



Alternatively, one can use an iterative solution of (29) diverging to infinity as 

^ = l'°'(-i(i-"')'"(-i(i-''')'°(i(i -"'))))• "^i- 

where the number of the nested terms — ^ (1 — 77^) In (— ^ (1 — 77^) ln(. . .)) can be 
increased or decreased due to the desired precision. Expansion (31) also follows 
straightforwardly from the asymptotic properties of the Lambert W function [33]. Of 
course, in reality the energy grows until the temperature T from (15) turns to zero. In 
terms of /? and A;, the leading, A-independent coefficient of (31) can be deduced directly 
from (23): 

E k or/and ^ ^ 00. (32) 

The corresponding order parameter tends then to the /3-independent constant: 



1 + 

X(x) = xl — ^ — , k or/and ^ 00. (33) 

Note that at A; = 1 and very large nonhnearities, (32) and (33) constitute a solution 
of the initial normalized GL equation (19) for the bulk superconductor when no any 
boundaries are present: 

E\k=i^^, /3^oo, (34a) 

Xix)\k=,-1, P^oc. (346) 

This means that the effective potential (14) is so large in this limit that the confinement 
of the strip can not alter it. Equations (34) also show that the density with A; = 1 is 
a special one what will be vividly demonstrated for the negative de Gennes distances, 
subsection 3.2. 

In the opposite limit of the zero r] one arrives at the equations defining energies 
E in the linear regime [25,26,34], as expected. Other analytical asymptotics of the 



Interplay between boundary conditions and cubicity of Ginzburg- Landau equation 



9 



above equations is the case of the Neumann boundary condition, 1/A = 0, when the 
ground-state energy is an identical zero for any f3 and k, as it follows directly from (27a): 

E(A^oo,/3,k) = 0. (35) 

It means that for the Neumann strip the critical temperature does not depend on the 
cubicity of the GL equation and on the density of the superconducting carriers and is 
always equal to the bulk critical temperature Tc. As a final remark of this section, let 
us note that if one chooses as a unit of length not the strip width d but the coherence 
length ^(T) — \JW' j {2mE), then the changes Vtt^ d in the first argument of the 
elliptic functions and ^/tt'^E — )> 1 in all other appearances in (27) of the parameter ( 
from (24) allows one to find the critical width dc as a function of the de Gennes distance 
A and factor ri [21]. However, compared to the normalization implemented by us above, 
the cubicity /3 in such scaling does disappear from the corresponding equations since 
the role of the primary parameter in this case is played not by k but directly by rj that 
emerges as a result of the integration of the coefficient- free GL equation [21] 

x"ix) + xix)-xHx) = 0, (36) 

where the tilde denotes that the order parameter is measured in units of /3/\a\ [21]. 

3. Results and discussion 

Here, the results of the calculations based on the theory developed in the previous 
section arc presented and analysed from the mathematical and physical points of view. 
Since the obtained dependencies are different for the opposite signs of the de Gennes 
distance, the cases of A > and A < are considered separately. 

3.1. Zero and positive extrapolation lengths 

Figure 1 shows energy as a function of the parameter r] and positive dc Gennes distance 
A. Similar to the linear configuration [26], for all nonzero rj the energy decreases with 
the growth of the extrapolation length until in the Neumann limit, 1/A = 0, one gets 
the critical temperature of the film equal to its bulk counterpart, as it was stated in the 
previous section, equation (35). In turn, the growth of rj pushes the energy upwards 
and at 77 — > 1 one sees the divergence described by (30) and/or (31) with its steepness 
being larger for the smaller A. 

Evolution of the order parameter is shown in figure 2 where the panel (a) exhibits 
xix) with varying the de Gennes distance at the fixed f] = 0.5. Similar dependencies 
are characteristic for the other factors ri too. For the Dirichlet case, A = 0, due to the 
suppression of the order parameter at the edges, the nucleation of the superconductivity 
takes place in the middle of the strip. Growing extrapolation length allows the Cooper 
pairs to approach the interfaces with the corresponding decrease of the slope of the order 
parameter in the transverse direction, and, at the large A, the function x{x) almost does 
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Figure 1. Lowest energy E as a. function of tlie coefficient 77 and positive de Gennes 
distance A. At 77 approaching unity, the energy tends to infinity according to (31). 

not depend on x what means a uniform nucleation of the superconductivity across the 
film. 

Panels (b) and (c) of figure 2 provide a comparative analysis of the influence of the 
factor rj on the order parameter at the two extreme values of the extrapolation length: 
small, A = 0, panel (b), and very large, A = 100, panel (c). In either case, the symmetric 
function x{^) unrestrictedly grows when rj approaches unity. However, the x-shape for 
the Dirichlet case is much more convex as compared to the de Gennes distances close 
to the Neumann one, 1/A ~ 0, when, as panel (c) demonstrates, the order parameter 
stays flat with its magnitude close to zero for almost all values of 77 and starts to grow 
upwards only in the very neighbourhood of r/ = 1. The larger the de Gennes distance 
is, the sharper the rise is at the right edge of the interval from (28). 

According to (23), the index rj accommodates two variables (5 and k. To illustrate 
their separate influence on the temperature T, we plot in panel (a) of flgure 3 the 
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Figure 2. Order parameter x[x) [normalized to (2//?)^/^] as a function of (a) the de 
Gennes length A for = 0.5, and of the parameter rj for (b) A = and (c) A = 100. 
Note different vertical scale in panel (a) as compared to panels (b) and (c). 
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energy E in terms of these two factors for the extrapolation lengths A = and A = 1. 
It is seen that, as discussed in the previous section, for the zero value of either ^ or k 
the temperature does not depend on the second parameter. The energy increases and, 
accordingly, the temperature T decreases with the growth of f3 or/and k for any finite 
de Gennes distance. As stated in the Introduction, the temperature dependence on the 
cubicity is explained by the positiveness of the effective potential Veff{x), equation (14). 
In other words, the role of the nonlinearity increases with the temperature decreasing. 
The energy growth with the factor k is due to the fact that the decreasing temperature 
T leads to the increase of the density of the Cooper pairs that, as it follows from (2), 
(22) and (23), is determined, among other parameters, by k. As the calculations show, 
the rate of the energy growth along the positive k direction exceeds the one along the 
/3 axis: 

dE BE 

Ik > w '''' 

These two speeds of the energy change grow bigger with k and /3. 

As stated in the previous paragraph, the cubicity plays more prominent role when 
the temperature decreases. A natural question arises: when does this role become so 
dominant that the linear theory could not be applied and its predictions are not close 
anymore to the results obtained from the cubic GL equation? Practical importance 
of this problem lies in the fact that obtaining the solution of the linear equation is, 
generally, much easier than its cubic counterpart. To answer this question, in panels 
(b)-(d) of figure 3 the energy is plotted in terms of the extrapolation length A and 
cubicity /3 at several parameters k. It is seen that for the small densities, say, k = 0.1 in 
panel (b), the energy only weakly depends on the cubicity and, accordingly, the linear 
theory serves as a reliable substitution of the cubic GL equation in a wide range of the 
parameter (3. This is especially true for the Dirichlet case, since the energy change with 
the cubicity decreases for the smaller de Gennes distances: 

— E(Ai,/3,^)>— S(A2,/3,^), Ai>A2>0. (38) 

Growing density rig leads to the increase of the slope of the E-/3 characteristics, as a 
transition from panel (b) to (c) and, subsequently, to (d) demonstrates. Equation (38) 
holds true but the energy difference for the different extrapolation lengths gets smaller: 

-E{A,f3,ki) < ——E{A,(3,k2), ki>k2. (39) 



As a result, the interval of change of the coefficient f3 where the solution of the cubic GL 
equation can be replaced by its linear counterpart without noticeable loss of precision, 
shrinks and becomes less dependent on the de Gennes distance with growing density 
of the superconducting carriers n^. In fact, this conclusion is not unexpected since, as 
it was already stated above, the density is proportional to the square of the order 
parameter, equation (2), that, in turn, is defined by the factor k. At the same time, the 
same expression x^(x) enters the effective potential Veff{x), equation (14), what means 
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Figure 3. a) Lowest energy E as a. function of cubicity (3 and coefficient k for the 
Dirichlet case (upper surface) and for A = 1 (lower surface), (b)-(d) Energy i? as a 
function of the cubicity (3 and positive extrapolation length A for several values of the 
parameter k shown in each of the panels. 



that an increase of the superconducting density leads to the stronger deviation from the 
hnear regime. Let us note also that (39) leads to 

a^^(A.A*=)<0. (40) 

what is a consequence of the negative derivative with respect to the de Gennes distance 
and positiveness of the partial derivatives with respect to the cubicity and density, as 
it is exemplified in figure 3. Let us point out also that the approach to the limit from 
(32) for the large (3 is also clearly seen from panels (b)-(d) of figure 3. 
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3.2. Negative de Gennes lengths 

Before discussing the interaction of the negative de Gennes distance and cubicity, one 
needs to recall the linear hmit, ^ = 0, of the same geometry [17, 24-26, 34]. In this case, 
the ground-state energy for A < always lies below zero. The adjacent upper state 
possesses negative energy too but only if the condition 

A+ + A_ + 1< (41) 

is satisfied [26]. Absolute values of the negative energies increase with dwindling |A-|-|, 
and at the vanishing and, in general, different lengths A± they diverge as [11]: 

^± = -^, A+^-0, A_^-0. (42) 



TT 



± 



For the equal, equation (21), small negative de Gennes distances one has two almost 
degenerate states with their symmetric and antisymmetric order parameters localized 
mainly at the interfaces. For the smaller |A| the attraction of the interfaces gets 
stronger [26, 34], and at A — )■ —0 each of these two order parameters transforms into the 
corresponding superposition of the 5-functions with their origins at the interfaces [34] : 

X ev (x) ~ S{x - 1/2) ± S{x + 1/2), A ^ -0. (43) 

odd 

Thus, the negative extrapolation length transforms the interface into the attractive 
center what results in lowering the Cooper pair energy and its localization at the 
edge. The smaller the magnitude of the negative de Gennes distance is, the larger the 
attraction and the stronger the localization of the surface state are. Unlimited growth 
of the magnitude of the negative energy from (42) is extremely interesting from the 
practical point of view since it means, according to (15), an unrestricted growth of the 
temperature T. In the model of equations (9) and (10), the limit of the infinitely small 
negative Atot can be achieved, for example, by applying to the zero-field Neumann edge 
an appropriately directed weak electric field S. We also mention that (42) is a particular 
case of a more general property of the vanishingly small negative extrapolation length 
in any number of dimensions intensively studied recently by mathematicians [35-39] . 

As a first step in the analysis of the combined influence of the negative de Gennes 
distance A and cubicity /3 on the temperature T, let us note that, as it directly follows 
from the transformation properties of the elliptic functions with complex argument 
and/or parameter [27, 28]§, equations (22) and (27) remain valid for the negative energies 
E too. Alternatively, one can directly use (19) where for < it is necessary to put 
E = — l^"]. After some algebra one gets the order parameter in the form 



X- (x) = ^ ?7_C- nd (^C-X, ^Jl + r)'ij , 



(44) 



§ There is a typo in the Enghsh translation of Rcf. [27]; namely, the numerator of the right-hand side 
of equation 8.153.8 instead of the Jacobi elliptic sine should contain an elliptic cosine, as is the case in 
the Russian original. 
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Figure 4. Energy E a,s a function of the index ?/_ and negative extrapolation length 
A. 



where the factors ?7_ and are given as 
ri.=ri.{E,(3,k) = 



(45) 



(46) 



with their subscripts indicating that they correspond to the negative energies. 
Eigenvalue equation (27a) turns to 

C- (1 + 77!) cn(C/2, r/_) sn(C-/2, r]__) + j dn(C-/2, r]_) = 0. (47) 

Note that in this case the parameter rj^ strictly varies between zero and unity: 

< < 1. (48) 

Figure 4 depicts energy E as a function of the index t]_ and negative de Gennes 
length A. In the linear regime, t]- = 0, a peculiarity from (42) is clearly seen at the 
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vanishing A. Rising factor rj^ leads to the corresponding increase of the energy and in 
this way withstands the temperature growth produced by the negative A. Its pace of 
change becomes faster for the smaller |A|: 

At ?7_ = 1 the energy turns to zero what means, according to (15), a complete 
compensation of the influence of the negative length, 

T\r,_=i-n- (50) 

This equation holds for all extrapolation lengths. However, even though the de Gennes 
distance does not enter (50) explicitly, the temperature dependence on A is crucial here. 
To prove this, we state that 7/_ = 1 is one of the points on the semi infinite strip of the 
unit width {0 < ?7_ < 1, — oo < A < 0} in the neighbourhood of which the energy E as 
a function of the index rj- and the length A is represented in the separable form: 

£;(77_,A) = ^(,7!-1)7(A), rj^^l, (51) 

where, according to (47), the value of 7 is determined from equation 

/(7,A) = (52) 

with 

f{y,A) = 2yV2 cnQyV2^2V2^ snQ|/V^2V2) + 1 dn Q y V2 , 2 . (53) 

Then, from (44) and (51), the order parameter Xo{x) for T = Tc is written as 

Xo(:r) = y|y/^nd(y/^x,2V2). (54) 

Any other values of 7(A) from (51) that do not satisfy (52) are unphysical ones. To 
cxcmphfy this, in figure 5 we plot the order parameter, equation (54), for several values 
of 7 with the solid line representing the case when (52) is satisfied. Similar to the linear 
regime discussed above, rj- = 0, the probability of finding the Cooper pairs at the edges 
is larger than in the middle of the film. For the nonzero cubicity, rj- > 0, the effective 
potential from (14) pushes the energy upwards and only for the value of 7 satisfying 
(52) it exactly compensates the attractive influence of the interfaces with the negative 
A. The potential Veff{x) for any other 7 is cither too shallow to promote the energy to 
zero (see the dotted line in figure 5) or too strong (dashed line) with the corresponding 
temperature T lying already below its bulk critical counterpart what in either case 
means that the temperature T(?7_,A) in the vicinity of = 1 is not a continuous 
function what, obviously, eliminates these two unphysical 7 from the consideration. 

Similar to the case of the positive extrapolation lengths, the coefficients /3 and k 
are compacted into the single factor 77- , according to (45). In terms of /3 and k there are 
three distinct situations: A; < 1, A; = 1 and k > 1. Prom the discussion above it directly 
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Figure 5. Order paranieter Xo{x) (in units of ^^2/(3) from (54) for A = —0.5 and 
several values of the parameter 7 where the solid line is for 7 — 1.8341 . . . being a 
solution of (52), the dotted line is for 7 = 1, and the dashed line is for 7 = 3. Thin 
vertical line denotes a middle of the film. 



follows that the energy turns to zero for A; = 1 only with the additional demand on the 
cubicity to satisfy the relation 

/(f,A)=0. (55) 

From this equation a critical extrapolation length Aq at which the zero energy is reached 
[11] can be defined: 



Ao (/3) = 




It is a quantitative measure of the exact cancelling by the nonlinear GL term of the 
temperature increase produced by the negative extrapolation length and, as such, it 
allows to understand the interaction of the counterbalancing influences of the cubicity 
and the surface effects. Its limit for the small nonlinearities reads: 

Ao(/3) = -|, (57) 

This result is understandable since in the linear regime the energy for the large negative 
de Gennes distances lies only slightly below zero [26, 34] and, accordingly, small effective 
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potential (14) is sufficient to turn it to zero and subsequently to push into the positive 
territory. The magnitude of the negative Aq grows together with the cubicity until it 
reaches zero at 

Po = (2-'/' ) = ^^^^ = 13.7504, (58) 

4:71 

where r(a;) is F-function [27,29]. Vanishing Aq (/3o) means that all critical lengths 
with 13 > /3o, even though mathematically correct, are unphysical ones. In other 
words, the effective potential (14) with /3 — /3o (and, of course, k — 1) completely 
compensates an infinitely large temperature increase produced by the infinitely small 
negative extrapolation length. Any superconducting state with /3 > /3o and k — 1 has 
its temperature always smaller than the critical temperature of the bulk material, 

Ty>^o < Tc , (59) 

regardless of the surface contribution. Prom this point of view, the cubicity (3 of the 
GL equation as compared to the de Gennes boundary is a stronger factor affecting 
the properties of the superconducting film; however, as it will be shown in the next 
paragraph, a substantial role in this drastic energy increase is played by the de Gennes 
distance itself. Critical extrapolation length Aq is shown in panel (a) of figure 6 where 
the states above (below) the curve correspond to the temperature higher (lower) than 
T 

c- 

It might look surprising that the finite cubicity is able to bring into the area T < Tc 
the state with any vanishingly small negative A when, at /3 = 0, its energy tends to 
the negative infinity, according to (42). For example, a comparison with the other 
configurations shows that the critical extrapolation length Ao(-B), which is calculated 
from the linear GL equation for the quantum dot or antidot in the uniform magnetic field 
B, asymptotically tends to zero from below only for the intensities B approaching infinity 
[11]. To resolve this seeming paradox, one needs to analyse the form of the effective 
potential (14) for the small negative de Gennes distance that, in the linear case, (3 = 0, 
attracts the order parameter to the edges and, in this way, delivers the larger ability for 
the potential to repulse the Cooper pairs from the walls. Thus, the extrapolation length 
tending to zero from the left produces at nonzero cubicity two contradicting to each 
other effects: on the one hand, it pushes the energy downwards by localizing the order 
parameter at the interfaces and, on the other hand, as a consequence of the previous 
action, it simultaneously provides the effective potential Veff{x), equation (14), with 
the larger ability to repel the function x(^) from the film edges with the corresponding 
increase of the energy. The smaller the magnitude of the negative de Gennes distance 
is, the larger the negative energy becomes and the stronger localization at the surfaces 
is but, at the same time, just due to this stronger localization, the faster the effective 
potential Veff{x) grows with the small increasing (3. In other words, potential (14) in 
its effort to push the energy upwards and the order parameter away from the edges 
is actually aided by the small negative length A that facihtates this energy growth 
by pulling the order parameter to the interfaces with the corresponding increase of its 
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Figure 6. (a) Critical extrapolation length Aq, equation (56), and (b) corresponding 
to it zero-energy order parameter from (60) versus the cubicity /?. 
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amplitude and in this way subjecting the function x{^) to the opposite influence of the 
effective potential. This fine interaction between the cubicity and de Gennes distance 
results ultimately in the drawing out of the energy from any arbitrary large negative 
values. 

Concerning equation (56) and figure 6(a), it is necessary to elaborate the following 
remark already mentioned in passing while discussing (58). Since the Jacobi elliptic 
functions - in particular, sn(M, k), cn{u, k) and dn(M, /c), which enter (56), - are doubly 
periodic functions of its argument for each Aq from figure 6(a) there exists a countably 
infinite set of its counterparts with the same value achieved at the larger cubicities j3. 
However, physically all these critical extrapolation lengths for (3 > f3o are spurious ones. 
This can be seen, for example, from the fact that their corresponding order parameters 
exhibit infinite discontinuities at some points inside the film, —l/2<x< 1/2. To save 
space, we do not plot these unphysical functions here presenting instead in panel (b) of 
figure 6 an evolution of the zero-energy order parameter 

Xo(^)=nd|./?a;,2V2) (60) 




2 



with the growing cubicity. It exhibits a smooth transformation from the flat Neumann 
shape at /3 = 0, when, according to (56) and (57), the critical extrapolation length is 
equal to the negative infinity, to the linear superposition of the two 5-functions, equation 
(43), at /3 = /3o when the critical de Gennes distance tends to zero from the left. For any 
cubicity larger than the critical one, (3 > (3o, these infinities migrate from the edges into 
the middle of the film turning in this way the corresponding solutions into the spurious 
ones. 

Next, consider the cases of A; 7^ 1. When the integration factor k is less than unity, 
k < 1, the corresponding solution always exists with its energy lying below zero for all 
cubicities. For the large nonlinear term, /3 — )■ 00, it ceases to depend on the cubicity 
approaching asymptotically a solution of the equation 



(1 + M o Ho\ T^2^ ^/i^ 



l-k^^ ' I 2 V 1 - A;2 ' / I 2 V 1 - A;2 



+ ^dn|^^y^^,^rTFj =0, 0</c<l, /3^oo, (61) 

that is derived from (45) - (47). For the smaller densities the saturation energy 
(temperature T) is smaller (larger). Corresponding order parameter is given then as: 



27r2|E|A;2 



^^"^^ ^ V /3(1 - P) ^ ^) ' 0<A;<1, ^^00, (62) 

what means that the effective potential Ve//(x) in this limit becomes /3 independent: 



VefAx) ^ ^^y]f nd^( y X, + ¥ I , 0<A;<1, /3^oo, (63) 
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and it is too weak to subdue the state sufficiently enougli for its temperature to lower 
to Tg. On the other hand, for k > 1 the energy can not grow more than the threshold 
value Eth defined by the cubicity and density: 



as it directly follows from (45) and (46). In fact, the state ceases to exist at the energies 
smaller than Eth since (47) in this case contains the Jacobi functions of the form 



that, similar to their trigonometric counterparts [40], possess infnitely many oscillations 
for E — >■ Eth- Upon its growth with the cubicity, the eigenvalue of (47) approaches one 
of the zeros of these oscillations and merges with it after which the two corresponding 
rea/ energies transform into the two comj^/ex conjugate values. Since the present research 
deals with the real energies and order parameters only, this merger is considered 
as a disappearance of the corresponding state. Physically, a destruction of the 
superconductivity for these parameters is explained by the fact that the decreasing 
temperature T can not support the fixed large density rig of the superconducting carriers 
with its concentration near the film interfaces. 

To exemphfy the results obtained above, figure 7 depicts energy as a function of the 
cubicity for several miscellaneous combinations of the negative extrapolation lengths and 
densities. It is seen that for k < 1 (dashed line) the energy saturates with the increasing 
f3 to its asymptotic value calculated from (61) while for k > 1 (dash-dotted curve) the 
state ceases to exist on its energy approach to Eth- A comparison of the curves with 
k = 1 and different magnitudes of the negative de Gennes distances (solid, dotted and 
two dash-dot-dotted lines) vividly manifests a convergence of the zero-energy cubicity 
to its critical value from (58); namely, for the large enough negative A only the small 
cubicity is needed to bring the energy to zero (dotted line in figure 7) while for the very 
small magnitudes of the de Gennes distance, for example, depicted by the two dash- 
dot-dotted lines, the energies of the corresponding states that differ by four orders of 
magnitudes in the hear regime, P = 0, cross zero at about the same /3 being very close 
to /So- After passing zero, the energy continues to rise with its order parameter given by 



As it follows from (34 a), for the very large nonlinearities, it diverges to the infinity with 
P/tt^ being its principal term, and the order parameter tends to unity, equation (346). 
Thus, since the concentration with A; = 1 for the large cubicities corresponds to the 
boundary-free situation, the state with only this density is allowed, for the negative de 
Gennes distances, to cross the bulk critical temperature with the increasing /3. Evolution 
with the cubicity of the order parameter x{^) corresponding to the solid line of figure 7 
is shown in figure 8. It is seen that its upward concavity decreases with the growing 



E < Eth = 



k>l, 



(64) 





(65) 
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Figure 7. Energy as a function of tlie cubicity /? for several values of the density 
and negative extrapolation length where the dotted curve is for A = — 2 and k — 1, 
the solid line is for A = —0.2 and k ~ 1, the dashed line is for A = —0.2 and k = 0.8, 
the dash-dotted line is for A = —0.2 and k = 1.2, and two dash-dot-dotted lines are 
for A: = 1 and A = —0.01 (upper curve) or A —0.0001 (lower curve). Thin horizontal 
line denotes zero energy. 



parameter /3 until, at the very large nonlinearities, it becomes, according to (346), 
completely flat. 

Finally, we provide a very instructive comparison of the above results for the 
temperature T passing through the bulk critical temperature with the predictions 
derived from (19) where the energy E is explicitly set to zero: 

X" (x) - /3x' (x) = 0. (66) 
Its corresponding eigenfunction x^*^^ (x) reads: 



AO), 



//3 



r"'[x) = knd{ k\l'-^x,2^/^] . (67) 
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0.50 



Figure 8. Order parameter as a function of tlie cubicity (3 for A = —0.2 and 
k = 1. Corresponding energy is depicted by the solid line in figure 7. 



Then, it directly follows from the boundary condition 



f 



.A 







that the corresponding critical extrapolation length A^^^ is: 




{61 



(69) 



A comparison with the previous derivation shows the difference between the two results; 
namely, while (56) is obtained under the strict condition that the integration factor k 
is unity and cubicity /3 obeys (55), the claim of (69) is that the zero energy can be 
achieved at any arbitrary values of the positive parameters /3 and k provided that their 
product satisfies (68). It is easy to understand the origin of this wrong conclusion; 
namely, when the energy is bluntly put to zero in (19), resulting equation (66) and 
its solution (67) know nothing about the quasi-peculiarity in (46) and its elimination 
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by the ansatz from (51). Equations (44) - (47), (55), (56) and (60) for the negative 
E contain much more information about the temperature behaviour close to than 
their zero-energy counterparts (67) - (69). Due to their hmited knowledge, the latter 
equations erroneously advice us that the zero energy can be achieved at any magnitudes 
of (3 and k if the product k'^(3 obeys (68) while in the reality the film temperature is 
equal to the bulk critical temperature T^. for the unity factor k only with the cubicity (3 
satisfying (55). 

4. Concluding remarks 

Exact solutions of the GL equation for the superconducting film with its edges 
subjected to the Robin-type boundary condition (12) revealed strong dependence of 
the temperature T on the interaction between the de Gennes distance and cubicity 
with, in addition, very substantial influence of the superconducting carriers density n^. 
Depending on these parameters, a rich spectrum of the temperature T is obtained and 
analysed in detail. For example, the interplay of the negative extrapolation length and 
nonlinearity causes the decreasing with the growing cubicity temperature T to cross 
its bulk critical value at the unique carrier concentration only while for the smaller 
densities its stays above for all cubicites, and the superconductivity is destroyed by 
the increasing coefficient /3 for the larger Gooper pairs concentration. The uniqueness 
of this density stems from the fact that only for it the increasing cubicity transforms 
the solution of the GL equation for the film into the one corresponding to the bulk 
superconductor. 

There are several possible ways of varying the de Gennes distance and magnitude of 
the nonlinear GL term. Experiment suggests that the value of the extrapolation length 
can be controlled by the sample preparation procedure [15]. Even more promising 
from this point of view is its manipulation by the applied gate voltage, according to 
(9) and (10). Miscellaneous aspects of the influence of the electrostatic flelds on the 
superconducting fllms have been the subject of the experimental investigations for more 
than flfty years [41-47]; in particular, the change of the critical temperature under 
nonzero £ has been observed [41,45,47]. However, despite its long history of research, 
electric fleld effects on superconductors are still far from being understood [48-53]. 
With regards to the cubicity ^, the methods of its variation follow straightforwardly 
from deflnition (4) for the pure metals [7, 8] and, especially, its alloy modiflcation [6, 7, 9] ; 
namely, different impurity doping drastically changes the mean free path I that becomes 
much smaller than the coherence length ^(0) thus leading to the substantial increase 
of the nonlinear influence. In addition, the range of validity of the GL theory is much 
wider for the 'dirty' superconductors [7]. Accordingly, it is believed that applying the 
gate voltage to the fllms composed with the different concentration of impurities, one 
can experimentally observe the predicted phenomena. 
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